0254-6124/2015/35(5)-525-08 Chin. J. Space Sci.。 空间 科学 学 报 


Zhang Xiangyu, Zhao Hua. Simulation of induced step-transient magnetic field with a lunar homogeneous model (in Chinese). 
Chin. J. Space Sci., 2015, 35(5): 525-532, doi:10.11728/cjss2015.05.525 


月 球 均匀 模型 阶 跃 感应 磁场 模拟 计算 
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摘 要 针对 月 球 均匀 模型 , 利用 电磁 感应 理论 , 对 行星 际 磁场 阶 跃 扰动 产生 的 感应 磁场 进行 了 模拟 计算 . 设 磁 导 
率 恒 为 jo, 对 于 一 些 特定 的 电导 率 值 , 给 出 了 月 表 磁 场 分 量 和 总 场 从 路 变 到 重新 达到 稳 态 的 变化 过 程 ; 在 一 条 经 线 
上 的 赤道 附近 、 中 纬 地 区 和 极 区 各 选取 一 个 测 点 , 利用 Laplace 北 变 换 的 数值 公式 给 出 了 磁场 瞬 变 响应 函数 和 磁 
场 分 量 及 总 场 在 不 同 电 导 率 情况 下 的 变化 趋势 . 计算 结果 表明 , 此 研究 方法 可 行 , 结果 合理 . 经 过 足够 长 的 时 间 , 与 
外 磁场 牙 变 方向 垂直 的 磁场 分 量 将 会 消失 , 而 与 之 平行 的 磁场 分 量 将 与 外 场 趋 于 一 致 . 在 外 场 发 生路 变 的 时 刻 , 未 
道 附近 能 够 测 到 的 最 大 磁场 分 量 值 约 为 13.65 nT, 极 区 附近 能 测 到 约 2.71 nT; 在 中 低 纬度 和 极 区 , 平行 分 量变 化 
的 过 程 显著 不 同 , 反映 出 实际 探测 中 , 选取 在 月 表 不 同位 置 的 磁 强 计 将 可 能 监测 到 完全 不 同 的 磁场 变化 曲线 形态 . 
关键 词 正 演 , 球 模型 , 电磁 感应 , Laplace 逆 变 换 
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Abstract With a homogeneous model of the Moon, the induced magnetic field caused by step- 
transient disturbation of interplanetary magnetic field is simulated by the electromagnetic induction 
theory. As the range of conductivity assigned, the variation from step to calm of both the vectors and 
magnitudes can be obtained with a fixed permeability /0. Then the magnetic characteristic transient 
response function, the field components and total field at 3 different points selected in a longitude 
line can all be given in different conductivities by numerical inverse Laplace transformation. The 


calculation provides a reasonable result that components perpendicular to the external step-transient 
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field will disappear in a long-enough time, whereas the parallel component will become unity with it. 


When the external field step varies, the maximum field value which can be detected in the vicinity 


of the equator and the polar region are 13.65nT and 2.71 nT, respectively. The processes of parallel 


component change are remarkably different in low latitude and polar region. 


This phenomenon 


reflects that magnetometers in different sites of lunar surface may lead to different field curves. 


Key words Forward modeling, Spherical model, Electromagnetic induction, Inverse Laplace 


transformation 


0 引言 


月 球 磁 场 是 当前 空间 探测 的 研究 热点 , 而 行星 际 
空间 的 扰动 磁场 是 很 好 的 研究 工具 .月 球 没有 全 球 
性 的 偶 极 磁场 , 因此 行星 际 磁场 一 旦 发 生变 化 , 再 经 
过 足够 长 的 时 间 , 就 可 以 完全 穿 透 月 球 内 部 . 这 一 过 
程 中 , 在 月 表 监 测 到 的 磁场 大 小 和 方向 变化 数据 携带 
有 月 球 的 分 层 结构 、 各 层 厚 度 、 磁 导 率 、 电 导 率 、 横 
向 或 纵向 可 能 存在 的 非 均 匀 性 变化 以 及 由 此 反映 出 
的 可 能 的 岩石 分 布 和 组 成 成 分 、 地 质 特点 、 物 理化 
学 性 质 等 一 系列 深部 信息 .因此 通过 在 月 表 测量 的 
相关 数据 , 再 利用 某 些 方法 和 手段 即 可 对 这 些 信息 进 
行 深入 分 析 与 研究 . 

1969 年 美国 Apollo 12 和 Explorer 35 磁 强 计 
分 别 在 月 面 和 环 月 轨道 对 几 个 行星 际 磁 场 阶 路 扰动 
事件 进行 观测 ， 同时 记录 了 磁场 三 分 量 的 变化 ， 这 
是 首次 在 月 表 对 磁场 进行 的 直接 探测 , 随后 Apollo 
14，15，16 和 俄罗斯 Lunokhod 2 磁 强 计 也 对 月 表 
进行 了 磁场 测量 .探测 结果 表明 , 月 球 可 能 存在 两 
种 结构 .一 种 为 两 层 结构 ， 内 层 半径 为 0.95 Rn (月 
球 半 径 ), 外 部 壳 层 电导 率 极 低 , 几乎 为 零 , 可 忽略 不 
计 , 外 磁场 知 发 生变 化 瞬间 即 可 穿 透 月 壳 ; 用 两 层 模 
型 与 取得 的 数据 反 演 结果 进行 最 小 二 乘法 理论 对 比 ， 
得 到 的 最 佳 电导 率 匹 配 值 为 1.5x10-“Sm- 1， 另 一 
种 为 近似 三 层 结构 ， 由 内 向 外 月 核 、 月 慢 、 月 壳 的 
厚度 分 别 为 0.6 Rw，0.37 Rm,，0.03 Rm; 电导 率 分 别 
为 10-?9，10-4 和 10-2Sm-1; 部 分 研究 得 到 的 一 
维 谢 面 不 完全 一 样 , 但 基本 上 电导 率 的 变化 范围 都 
在 10-7~10-5S.m-l 到 10-2 ~ 1S.mr-l2-4; 通过 
数据 分 析 得 到 月 球 平均 磁 导 率 为 (1.03 士 0.13)Hpo,， 如 
果 内 核 居 里 温度 高 于 800°C, 则 推断 得 到 的 核 外 磁 导 


算 设 月 球 磁 导 率 恒 为 yo; 但 对 于 目前 推断 出 的 这 两 
种 分 层 结 构 而 言 , 各 层 电导 率 可 能 存在 几 个 数量 级 的 
差别 . 由 于 本 文 只 研究 均匀 模型 , 所 以 将 电导 率 范围 
设 定 为 10-” ~ 10“Sm- ,接近 中 间 层 的 值 ; 行星 
际 磁场 的 变化 幅 值 一 般 不 超过 10nT, 因此 本 文 计算 
设 定 外 磁场 在 上 = 0 时 有 10nT 的 路 变 . 

Apollo 系列 磁 强 计 均 只 在 月 表 进 行 了 单 点 探测 ， 
未 来 在 技术 手段 可 靠 的 前 提 下 可 以 在 月 表 进行 多 点 
联 测 , 本 文 即 是 针对 不 同 均一 模型 , 通过 计算 给 出 部 
分 测 点 的 磁场 变化 值 ， 更 进一步 的 目标 应 该 是 能 够 
通过 月 表 不 同 测 区 实测 数据 的 对 比 和 联合 反 演 得 到 
更 高 精度 的 月 球 模型 , 而 本 文 的 正 演 计算 可 为 反 演 提 
供 一 些 前 期 的 理论 验证 . 


1 原理 及 研究 方法 


月 球 与 地 球 和 太阳 的 相对 位 置 是 随时 间 而 变化 
的 .为 简单 起 见 , 本 文 忽略 这 种 变化 , 研究 全 空间 磁场 
跃 变 作 用 于 月 球 的 情况 , 通过 计算 得 出 月 表 磁 场 的 变 

根据 电磁 感应 理论 , 外 部 磁场 的 变化 可 以 在 球体 
内 部 激发 涡 电 流 , 从 而 产生 二 次 感应 磁场 , 其 幅 值 和 
相位 与 外 源 磁场 相 比 会 发 生 一 定 偏差 . 实际 探测 过 
程 中 , 在 球体 表面 上 测 到 的 值 是 外 源 磁场 和 感应 磁场 
的 县 加 值 ; 由 于 二 次 场 会 随 着 时 间 的 推移 逐渐 衰减 直 
至 最 终 消 失 , 所 以 经 过 足够 长 的 时 间 之 后 , 球体 表面 
得 到 的 磁场 最 终 会 与 外 源 场 趋 于 一 致 四. 
1.1 ”模型 假设 

采用 如 下 模型 假设 , 即 月 球 是 一 个 各 向 同性 的 不 
分 层 均匀 介质 球体 , 不 考虑 电导 率 的 横向 变化 和 不 同 
区 域 地 形 高 差 变 化 ; 月 球 具 有 单一 电导 率 c 和 单一 磁 


率 为 (1.05 上 0.14)700,46-7， 由 此 可 见 , 与 真空 磁 导 
率 Jo 相 比 , 全 月 磁 导 率 没有 太 大 变化 , 故 文中 的 计 


导 率 /. 这 样 从 Maxwell 基本 电磁 方程 组 出 发 , 能 够 
推导 出 磁场 互 和 磁 矢 势 A 所 满足 的 方程 , 再 加 入 定 
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解 条 件 , 即 可 求解 该 方程 . 

利用 Laplace 变换 对 磁场 扩散 方程 进行 求解 , 并 
对 第 三 类 边界 条 件 进行 同样 处 理 ， 在 保证 精度 的 前 
提 下 , 通过 数值 方法 给 出 不 同 电 导 率 情况 下 月 表 磁场 
三 分 量 的 变化 过 程 , 进而 得 到 具体 部 分 测 点 的 分 量变 
化 , 并 进行 分 析 和 解释 . 

1.2 ”球体 电磁 感应 理论 

描述 球体 对 任意 外 磁场 输入 的 响应 ( 见 图 1), 定 
义 磁 矢 势 A, 其 满足 VxA 一 B 和 VvV.A=0. 
寻求 人 Bb(t) 这 样 一 个 阶 牙 变化 的 响应 ， 其 中 b(t) 
满足 下 列 条 件 : b(t) = 0 (t < 0) 且 5() 是 一 个 定 
值 (t > 0); 对 于 阶 路 磁场 ,有 b(t) = 1(t > 0). 由 于 控 
制 方程 是 线性 的 , 因此 一 般 的 输入 响应 亦 可 通过 从 加 
求 得 向 . 

AB。 方 向 设 定 为 球 坐标 系 (7,0,w) 的 轴线 ， 
即 Descartes 坐标 系 的 z 轴 方 向 ， 当 电导 率 呈 球 对 
称 时 , 瞬 变 磁场 响应 没有 2 分 量 , 因此 4 = 4ew 并 
且 对 所 有 分 量 均 有 0/9% = 0. 在 这 些 假定 条 件 下 ( 忽 
略 位 移 电 流 ), Faraday 电磁 感应 定律 、Ampere 定律 
和 Ohm 定律 共同 给 出 了 静 磁 场 中 的 Maxwell 电磁 
方程 组 及 本 构 关 系 , 消 项 后 可 以 得 到 如 下 无 源 磁 场 扩 
散 方程 : 

OA(r, 9,t) 
At 
其 中 , 4 和 o 分 别 为 球 内 介质 的 磁 导 率 和 电导 率 . 在 
外 加 激励 磁场 的 情况 下 , 介质 内 的 感应 磁 和 天 势 应当 满 

足 此 方程 . 

在 上 文 已 经 假定 在 全 空间 中 磁 导 率 满足 / = jo 
这 样 , 当 + > 0 时 , 磁场 在 球体 表面 边界 上 应 当 保 
持 连续 , 即 A 和 “4 在 + = ,处 连 续 球 心 满 
足 A(0,0,)=0; 球 内 满足 初始 条 件 A(7,0,0)=0. 在 


V*?A(r,0,t) = uo (1) 


ABeb(t) /nT 


0 2 5 10 20 


图 1 外 磁场 ABeb(t) 
Fig.1 Exterior magnetic field AB.eb(t) 
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球 外 由 于 是 真空 , 故 o = 0. 对 Laplace 方程 取 分 离 
变量 解 , 可 以 得 到 
A=AB。 (5) b(t) sin 0+ 


A f(t) sin0. (2) 


式 (2) 右边 第 一 项 由 给 定 的 外 部 磁场 变化 函数 b(t) 
决定 ; 第 二 项 即 所 求 的 外 场 瞬 态 响应 函数 , 当 7 一 oo 


和 上 一 co 时 将 会 消失 . 
将 7 归 一 化 后 对 公式 进行 整理 和 变换 , 可 以 得 到 
oA 3 5 
有 一 24 十 5IAB。 sin0 b(t)]. (3) 


由 于 磁场 在 7 = Rw 处 连续 , 式 (3) 即 是 对 球体 内 部 
磁场 进行 计算 时 的 边界 条 件 , 其 形式 为 第 三 类 混合 边 
界 条 件 . 令 

g(7,t) = A/ABo。 sin0 


并 用 G(r, s) 表示 g(7,t) 的 Laplace 变换 , 则 式 (1) 可 
以 写 为 


[总 69) ——G| = HasG (4) 
边界 条 件 式 (3) 可 以 转换 为 如 下 形式 : 
OG 3 
一 2G 十 30(5). (5) 


式 (5) 在 7 = Rw 处 成 立 , 其 中 , 5(s) 为 b(t) 的 Laplace 
变换 . 在 球 心 处 还 有 男 一 边界 条 件 , 即 


G=0. (6) 


这 样 , 在 给 定 电导 率 c 和 函数 b() 之 后 , 即 可 通 
过 数值 积分 获得 球体 任意 区 域 ( 含 球面 ) 的 G(" s); 
用 数值 方法 求解 方程 Laplace 逆 变 换 , 即 可 研究 系统 
的 特征 瞬 态 响应 函数 f(t) 及 其 性 质 . 求 出 磁 和 天 势 A 
后 , 将 其 用 B= V x 4 转化 为 磁场 值 即 可 . 


2 Laplace 道 变 换 


在 计算 过 程 中 给 定 二 维 函数 Laplace 变换 的 生 
成 函数 G(7, s), 需求 出 与 之 相对 应 的 决定 函数 g(7,). 
但 实际 上 大 部 分 函数 的 Laplace 逆 变 换 都 只 能 求 数 
值 解 包括 本 例 )， 在 Laplace 逆 变 换 的 定义 中 , 生 
成 函数 在 解析 半 平 面 一 致 收敛 ,而 积分 是 沿 着 这 个 
半 平 面 内 某 条 与 虚 轴 平行 的 直线 进行 的 ， 因此 对 于 
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其 数值 积分 不 能 完全 按照 常规 的 实数 积分 域 进 行 处 
理 , 而 需要 在 此 基础 上 进行 一 些 修 正 . 参照 Dubner 
和 Abatel9| 及 Durpinfol 的 方法 进行 如 下 变换 . 

对 于 本 文 的 生成 函数 , Laplace 道 变 换 的 定义 公 
式 如 下 : 


a 


27i 


Qa 十 icoe 
去 | 本 estG(7, s)ds. (7) 
其 中 , 复 变 量 s = a+iw, g(7,t) 的 Laplace 变换 在 半 平 
面 Re(s) > c 上 一 定 存 在 , 这 里 c 由 生成 函数 G(m s) 
决定 , 且 在 此 平面 内 G(7, s) 为 解析 函数 . 
将 式 (7) 变换 为 


9eat 


g(7,t) = | Retce 5)} coswt dw. (8) 


其 中 , a 为 大 于 c 的 任意 实数 . Dubner 和 Abatel] 对 
式 (8) 应 用 梯形 近似 公式 得 到 生成 函数 的 数值 反 演 
近似 公式 , 并 对 反 演 误差 进行 讨论 ，Durbinl10 还 对 
文献 [9] 的 算法 进行 了 一 些 改进 . 采用 文献 [9] 的 公 
式 并 用 在 式 (8) 的 二 维 函 数 中 , 得 到 如 下 反 演 公式 : 


2ert | 1 
= (ee a) 十 
元 


g(r,t) 3 


nle(nat 各 )] (多 )} + 


其 中 , i 为 虚数 单位 , 了 为 需要 计算 的 数据 t 的 长 度 ; 


bE = >》， e277T [g(r, 2nT + t+ 


= 


e2tg(r, 2nT — t)]. (10) 


其 中 , 参数 a > c, 7 > tmax. 式 (10) 中 一 旦 确定 了 了 ， 
利用 梯形 公式 中 的 步 长 r/ 开 可 以 近似 求 得 9 为, 这 
样 通过 选取 合适 的 a 代入 计算 , 可 以 使 得 误差 达到 预 
设 精度 . 为 了 确定 误差 的 大 小 , 需要 研究 当 t 较 大 时 
函数 g(7,t) 的 性 质 . 上 文 已 假定 生成 函数 在 一 致 收敛 
平面 Re(s) > c 内 无 奇 点 , 故 有 


g(r,t) < Ct™. 


其 中 C 为 常数 , m 为 任意 非 负 整数 . 
假定 对 任意 t 均 满 足 g(7,t) < C, 则 式 (10) 可 

以 改写 为 

a(t_T) cosh(at) 


及 < : 。 
oe sinh(a7 ) 


(11) 
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由 式 (11) 可 知 , 当 t 二 工时 , 误差 可 能 会 达到 不 
可 控 的 程度 , 因此 计算 的 有 效 范围 为 


0<t<T7/2. 


对 于 数值 公式 (11), 误差 限 为 Ce "7 选 定 了 后 再 
指定 a, 则 对 于 (0, 7/2) 范围 内 的 反 演 , 误差 可 以 控 
制 在 C x 4.54 x 10-5 以 内 . 


3 计算 结果 及 分 析 


本 计算 采取 的 模型 参数 设置 为 , 月 球 磁 导 率 1 = 
Ho = 4r x 10-7Hm-l; 月 球 电 导 率 o = 10-” ~ 
10-5S.m-1; 月 球 半径 Ru = 1738km; 计算 时 长 t = 
20s. 对 于 阶 跃 变 化 的 外 磁场 , 在 球面 (7 = 1) 上 , 频 
域内 满足 定 解 条 件 式 (4) 的 解 为 


G1(s) =G!(1, s) = 


[kcosh(k) — sinh(%k)]. 


12s(ek — ek) 


其 中 ,有 = Yhos5. 由 于 了 > 2tmax, 设 卫 = 50000， 
aT = 12, 则 a = 2.4 x 10-4, 0 < 5.1125 x 
x10-10. 式 (9) 截断 项 取 到 k = 300. 给 定 这 些 参 数 
能 够 保证 计算 达到 预定 精度 . 

图 2~5 分 别 给 出 了 球面 测量 磁场 的 三 分 量 及 
总 场 随 时 间 的 变化 趋势 ， 球 内 电导 率 均 为 o = 
10 “Sm 1. 可 以 看 出 , 在 外 磁场 发 生路 变 的 瞬间 , 球 
面 立 即 产生 响应 , 计算 得 到 最 大 的 B; 分 量 响 应 值 约 
为 5.44nT, 出 现在 Ozz 面 与 球体 相交 的 两 条 经 线 南 
北 中 纬度 地 区 ; 最 小 的 Bu 分 量 响应 值 约 为 -5.44nT， 
出 现在 与 两 条 经 线 对 称 的 南北 中 纬度 地 区 ; 最 大 By 
分 量 响应 值 和 最 小 分 量 响应 值 与 B。 相同, 且 出 现 
的 区 域 与 Bs 相应 区 域 的 纬度 相同 , 但 经 度 向 东 偏 移 
了 7/2; 随 着 时 间 的 推移 , 月 表 磁 场 与 月 内 地 质 体 相 
互 作用 , 导致 Bs 和 B 分 量 显著 衰减 , 最 大 响应 值 
到 11.4s 左右 已 降 至 约 0.24nT, 到 19s 时 更 是 仅 有 
不 到 0.04nT, 经 历 足 够 长 的 时 间 后 将 衰减 殖 尽 ， 赤 
道上 的 B; 和 By 分 量 始终 为 零 . 

由 于 假定 的 阶 牙 磁场 沿 轴线 z 的 方向 变化 , 所 
以 B; 分 量 的 变化 与 B; 和 By 显著 不 同 . 这 种 差异 
不 仅 体 现在 场 值 的 大 小 , 还 体现 在 分 量 在 球面 的 对 称 
性 上 . 从 图 4(a)~(f) 可 以 看 出 ,，B; 分 量 在 极 区 的 初 
值 仅 为 2.71nT, 而 在 赤道 附近 初 值 可 达 13.65nT, 随 
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(a) 


t=1.52s 


t=0.38s 


t=9.6s 
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t=15.2s 


图 2 o=10 Sm-! 时 月 表 Bs 分 量 的 变化 


Fig.2 B; at the lunar surface when o = 10 5S.m-! 


(a) (b) (0) 
上 
和 入 ， 
By/nT 
5.00 
3.75 


t=0.38s t=].528 


t=3.80s t=9.60s 


© -© Ol 


t=2.66s 2.50 


1.25 


t=15.20s 


图 3 o=10-5S.m-! 时 月 表 B, 分 量 的 变化 


Fig.3 By at the lunar surface when o = 10 ?SS.m-! 


着 时 间 推 移 , 极 区 附近 的 值 逐渐 增加 , 中 低 纬 处 的 值 
逐渐 减少 , 但 最 终 其 值 在 外 源 场 均 能 够 达到 10nT， 
且 B; 分 量 与 2 无 关 , 其 分 布 关于 赤道 南北 对 称 
图 2~5 都 是 在 电导 率 o = 10 Sm- 的 情况 
下 得 到 的 , 但 对 于 其 他 电导 率 , 球面 分 量变 化 的 总 趋 
势 是 一 致 的 . 由 于 磁场 在 岩石 中 的 衰减 速率 依赖 于 岩 
石 特性 , 因此 低 电导 率 条 件 下 磁场 衰减 快 , 高 电导 率 


条 件 下 磁场 衰减 慢 , 对 于 同样 的 外 源 场 , 高 电导 率 下 
测量 磁场 重新 达到 稳 态 所 需 时 间 将 会 比 低 电 导 率 下 
长 得 多 . 

这 种 变化 在 图 6~7 中 反映 得 相当 明显 . 在 球面 
上 选择 一 条 经 线 (p = 7/3) 上 的 位 于 北半球 的 三 
个 测 点 ,分 别 位 于 北半球 的 低 纬 (9 = 3r/7、 中 
纬 (0 = 2r/5) 和 高 纬 地 区 (0 = r/15). 图 6(a)~ (c) 
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t=1.52s 


t=3.80s t=9.60s 
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t=15.20s 


图 4 o=10-5S.m-! 时 月 表 Bz 分 量 的 变化 


Fig.4 B;, at the lunar surface when o = 10-5S.m-! 


t=3.80s 


t=9.60s 


t=15.20s 


图 5 o=10-5S.m-! 时 月 表 总 磁场 Bt 的 变化 
Fig.5 ‘Total field Bt at the lunar surface when o = 10-5S.m-1 


反映 了 各 个 测 点 的 三 分 量 及 总 场 的 变化 , 图 7(a)~(c) 
反映 了 轴 向 磁场 的 差 值 变化 . 可 以 看 到 , o = 10 ~ 
10“ Sm 时 , 三 分 量 仅 需 不 到 2s 甚至 1s 即 可 重 
新 达到 稳 态 , 但 是 当 o = 10 Sm 时, 这 个 时 段 将 
会 达到 8~10s 甚至 更 长 . B。 分 量 经 足够 长 时 间 后 将 
完全 等 同 于 外 磁场 , 但 极 区 与 中 低 纬 的 变化 趋势 也 不 


一 样 , 极 区 是 从 2~3nT 逐渐 增加 到 10nT, 而 中 低 纬 
和 赤道 是 在 t = 0 时 刻 先 突 增 至 某 个 峰值 , 然后 再 逐 
渐 减 小 至 10nT. 该 实验 结果 与 文献 [1] 中 Apollo 地 
面 实验 结果 是 一 致 的 , 说 明 本 文采 取 的 理论 算法 是 合 
理 的 . 但 是 当 o 较 小 时 ，B; 分 量 的 变化 经 过 一 定 振 
荡 才 逐渐 达到 外 场 值 , 表明 其 与 原 函 数 奇 点 有 关 , 这 
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图 6 球面 测量 磁场 三 分 量 及 总 场 (单位 nT 


随时 间 的 变化 趋势 
Fig.6 Three components and total field (Unit nT) 


We 


variations with time 


会 对 Laplace 数值 逆 变 换 的 公式 精度 产生 一 定 影 响 . 
为 了 压制 这 种 振荡 , 保证 磁场 平稳 变化 , 需要 控制 变 
换 公 式 和 误差 公式 的 参数 范围 并 进行 合理 选择 . 

图 8 描述 了 不 同 o 条 件 下 球面 控制 函数 g(t) 
及 外 场 特征 瞬 态 响应 函数 f(t) 的 变化 过 程 ， 根 
据 Laplace 逆 变 换 的 数值 反 演 公式 (9), g(t) 和 jb 
除 与 时 间 有 关外 , 还 与 球 半径 7 有 关 . 而 实际 的 测量 
不 可 能 深入 到 月 球 内 部 , 因此 在 球面 有 ” = 1, 两 个 
函数 均 不 受 测 点 位 置 的 经 纬度 影响 . 


假定 月 球 为 一 个 均匀 模型 , 在 球 外 空间 的 行星 际 
磁场 发 生 阶 路 变化 时 , 对 月 表 部 分 位 置 的 磁场 变化 进 
行 了 模拟 计算 , 得 到 如 下 结论 . 

(1) 月 表 磁 场 包含 行星 际 磁场 和 月 球 内 部 物性 参 
数 激发 的 变化 磁场 . 假定 磁 导 率 / 不 变 , 则 球 电导 率 
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一 G=105S.m 一 G=104S.m 一 aG=107Smr 


(a) 
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(b) 
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图 7 轴 向 > 磁场 差 值 随 时 间 的 变化 
Fig.7 Variation of the magnetic field in z-axis 


direction with time 


一 a=105S.m 一 oa=104Sm 一 a=107Smn 


图 8 不 同 电导 率 下 球面 控制 函数 g(t) 及 外 场 特 征 
瞬 态 响应 函数 f(t) 随时 间 的 变化 
Fig.8 Surface regulating function g(t) and 
characteristic transient response function f(t) 


with time under various conductivities 


0 越 小 , 变化 磁场 衰减 越 快 , 例如 o = 10- Sm- :时 
仅 需 2s 外 磁场 即 基 本 衰减 列 尽 ; 电导 率 越 大 , 变化 
磁场 衰减 越 慢 , 例如 o = 10-5S:m-l 时 需 约 20s 外 
磁场 才 基 本 衰减 殖 尽 ; 但 最 终 变化 磁场 都 会 消失 , 球 
面 各 点 的 总 磁场 均 将 无 限 趋 近 于 跃 变 后 的 外 部 磁场 ， 
这 符合 电磁 感应 定律 , 也 说 明了 Laplace 道 变换 数值 
反 演 方法 的 合理 性 和 有 效 性 . 
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(2) 对 于 球 对 称 模型 ， 南 北半球 的 磁场 变化 呈 
反对 称 情 形 ， 即 一 个 点 的 磁场 变化 和 经 度 与 其 相 
差 180°, 纬度 相同 但 方向 相反 的 点 的 磁场 变化 趋势 
相同 . 

(3) 对 于 均一 电导 率 模型 , 与 外 源 场 方向 垂直 的 
磁场 分 量 Bs 和 B, 最 终 都 将 消失 , 而 与 之 方向 相同 
的 分 量 B; 最 终 将 与 外 源 场 达到 一 致 . 但 是 平行 分 量 
在 中 低 纬 度 地 区 和 高 纬度 地 区 趋向 于 外 源 场 的 变化 
过 程 显著 不 同 , 因此 在 实际 测量 中 , 不 同 测 点 位 置 的 
选取 将 有 可 能 监测 到 完全 不 同 的 磁场 变化 形态 ， 这 
个 形态 差异 的 分 界 点 大 致 出 现在 南北 纬 38°. 

(4) 本 文 方程 求解 过 程 采 用 的 是 Laplace 道 变换 
数值 反 演 方法 . 已 有 文献 [9 一 10] 提 到 的 数值 算法 大 
多 研究 一 元 函数 , 由 于 用 到 了 复数 , 因此 由 精确 解 到 
近似 解 的 转化 需要 条 件 , 误差 限 的 控制 依赖 于 积分 函 
数 的 收敛 域 , 故 对 积分 函数 性 质 的 分 析 和 讨论 需要 并 
慎 对 待 ; 本 计算 中 的 方程 虽然 是 二 元 函数 , 但 在 计算 
过 程 中 仅 对 一 个 宗 量 上 进行 Laplace 变换 , 不 涉及 二 
维 情形 , 因此 对 一 元 函数 的 变换 过 程 完 全 适用 于 本 计 
算 , 对 计算 精度 和 收敛 条 件 的 讨论 亦 可 类 比 进行 . 

由 于 许多 研究 均 证 明月 球 存在 多 层 结构 , 因此 本 
文 计算 与 实际 探测 结果 可 能 会 存在 一 定 偏差 . 利用 
简单 均匀 对 称 模型 达到 预期 目标 , 证 明 布 设 月 表 磁 强 
计 阵 列 的 可 行 性 , 可 以 在 未 来 的 研究 中 改变 模型 结构 
和 参数 设置 , 例如 采用 分 层 、 偏 心 模 型 或 加 入 规则 、 
不 规则 结构 的 电磁 异常 体 , 考虑 地 层 结 构 的 横向 或 纵 
向 不 均匀 性 等 , 以 期 获得 更 接近 实际 , 具有 更 高 应 用 
价值 的 成 果 , 推进 不 同 测 点 差 值 与 垂 向 分 层 厚 度 之 间 
的 关系 及 反 演 异常 体位 置 和 尺度 等 方面 的 研究 , 为 月 
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球 探测 提供 算法 的 理论 依据 和 物理 过 程 分 析 . 
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